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SUMMARY 

A concept of system identification applied to high 
performance aircraft is introduced followed by a 
discussion on the identification methodology. Special 
emphasis is given to model postulation using time 
invariant and time dependent aerodynamic parameters, 
model structure determination and parameter estimation 
using ordinary least squares and mixed estimation 
methods. At the same time problems of data collinearity 
detection and its assessment are discussed. These parts of 
methodology are demonstrated in examples using flight 
data of the X-29A and X-31A aircraft. In the third 
example wind tunnel oscillatory data of the F-16XL 
model are used. A strong dependence of these data on 
frequency led to the development of models with 
unsteady aerodynamic terms in the form of indicial 
functions. The paper is completed by concluding remarks. 


NOMENCLATURE 

Only the main symbols are introduced. Other symbols are 
explained in the body of the paper. 


c 

mean aerodynamic chord, m 

b 

wing span 

c L 

lift coefficient 

c„ c m , C n 

rolling-, pitching- and yawing-moment 
coefficient 

k=o>yV 

reduced frequency 

t 

characteristics length, m 

p> q> r 

roll, pitch and yaw rate, rad/sec 

R 2 

coefficient of determination 

V 

airspeed, m/sec 

a, P 

angle of attack and sideslip angle, rad 


control stick deflection 

s, s, 

aileron and rudder deflection, rad 

5 C > 5 r S s 

canard, flap and strake deflection, rad 

V. 

thrust vectoring deflection in pitch and 
yaw, rad 

0 

unknown parameter 


1. INTRODUCTION 

The introduction of highly maneuverable and often 
inherently unstable aircraft has been presenting new 
challenges to aircraft identification and parameter 
estimation. These aircraft can perform rapid large 
amplitude maneuvers, often extended to the stall and 
poststall region where nonlinear and unsteady 
aerodynamic effects could be pronounced. This introduces 
a problem of determining how complex the model should 
be. Although a more complex model can be justified for 
more accurate description of airplane motion, it has not 
been clear in parameter estimation which relationship 
between model complexity and measurement information 
would be the best. If estimates for too many parameters 
are sought from a limited amount of data, a reduced 
accuracy can be expected or attempts to identify all 
parameters might fail. The high performance aircraft may 
also have more control surfaces moved through a flight 
control system than conventional aircraft. Such a system 
can introduce a close relationship between the deflections 
of various surfaces and at the same time can preclude 
maneuvers from being suitable for system identification. 
These characteristics can be reflected in an inability to 
estimate the effectiveness of individual control surfaces 
and to obtain accurate estimates of the remaining 
parameters. One of the reasons for these problems is 
related to the near linear relationship among several 
variables entering the model for various estimation 
techniques. This near linear relationship is often called 
data collinearity [1]. 

In recent years more attention has been given to the 
analysisof data obtained from dynamic wind tunnel tests. 
The frequency and amplitude dependency of oscillatory 
and ramp test data led to postulation of models with 
linear or nonlinear unsteady effects and subsequent 
parameter estimation in these models [2,3]. 

As follows from the above mentioned reasons and 

experience, a successful parameter estimation requiresthe 
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following. 


a) design of an experiment for obtaining data with high 
accuracy, high information content and low collinearity 
among measured inputs and outputs; 

b) determination of model structure which represents an 
adequate model for the aircraft under test in each 
maneuver analyzed; 

c) introduction of techniques which reduce the adverse 
effect of data collinearity on parameter estimates when 
severe data collinearity could not be avoided by 
experiment design. 

The purpose of this paper is to present a general approach 
to aircraft identification with emphasis on the above 
mentioned requirements for successful identification of a 
high performance aircraft model. The paper starts with an 
overview of system identification methodology followed 
by a discussion on data collinearity and biased estimation. 
In examples the flight data from experiments on X-29A 
and X-31A aircraft will be used. The forced oscillatory 
wind tunnel data on F-16XL model were selected for 
estimation of unsteady aerodynamic parameters. The 
paper is completed by concluding remarks. 

2. AIRCRAFT IDENTIFICATION 
METHODOLOGY 

When system identification is applied to an aircraft the 
equations governing its motion are postulated and an 
experiment is designed for obtaining time histories of the 
input and output variables. The equations of motion are 
formed by rigid-body force and moment equations 

mV +ma)xV =F q +F t + F(V,g>,u, 6) (i) 

I<j + <*) xlc*> = G(V,g>,u,0) 

and by a set of kinematic equations relating the Euler 
(attitude) angles and angular velocities. In eq. (1) m is the 
mass, I is the inertia matrix, V and to are the linear and 
angular velocity vectors, and u is the control vector. The 
vectors F c and F x represent the gravity and propulsion 
force, F and G aerodynamic force and moment 
respectively, and 0 is a vector of parameters which 
specify aerodynamic characteristics of the aircraft. 

For system identification, the aircraft state equations are 
completed by the output and measurement equations. The 
complete set of all these equations can be written as 

x=f[x(t),u(t),0], x(0)=x 0 
y =h[x(t),u(t),0] ^ 2> 

z(i) =y(i)+v(i), i = l,2,...,N 


where the state vector, x, is comprised of V, co and Euler 
angles, and u is the control vector of control surface 
deflections. The outputs, y, are the variables defining 
aircraft responses. The measured outputs, z(i), are 
corrupted by measurement noise, v(i), and the number of 
data points is N. 

Aircraft identification can be defined as follows: Aircraft 
identification is a determination, from input and output 
measurements, of a structure for F(0) and G(0) and 
estimationof unknown parameters, 0, in F(0) and G(0). 
In many practical applications, the structure of F (0) and G (0) 
is assumed to be known and aircraft identification is 
reduced to parameter estimation. A general approach to 
aircraft identification adopted at NASA Langley Research 
Center is shown in figure 1 in the form of a block 
diagram. Various steps in the procedure include model 
postulation, design of an experiment, data compatibility 
analysis, model structure determination and parameter 
estimation combined with collinearity diagnostic, and 
model validation. 

Model Postulation 

Model postulation is influenced by the type of selected 
maneuver intended for system identification and by a 
prior knowledge about aircraft aerodynamics. The 
aerodynamic forces and moments are expressed in the 
form of polynomials or polynomial splines as 

C.=C.(0) + E6x (3) 

j-1 

where C a is the aerodynamic coefficient, C a (0) is the 
value of the coefficient at initial steay-state conditions 
and the x. now represents input and output variables, 
their combinations and/or spline terms. The postulated 
model is then used in model structure determination and 
parameter estimation. For a model with unsteady 
aerodynamics the forces and moments can be formulated 
in terms of indicial functions [4,5] as 


c,(t) = C.(0) + f C, { (t-t;5M) T J- 5(T)dT ( 4 ) 


where £ is a vector of aircraft state and input variables 
upon which the coefficient C 4 depends, C^(t) is a vector 
of indicial functions whose elements are the responses in C a 
to unit steps in The indicial responses, C^, are 
functions of elapsed time (t-x) and are continuous 
single-valued functions of £(t). The indicial functions 
approach steady-state values with increasing values of the 
argument (t-x). If the indicial response is only a 
function of elapsed time, equations (4) is simplified as 
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C.(t) =0,(0) + jcjt-z) T j- 5(x)dT (5) 


When analytical forms of indicial functions are specified, 
the aerodynamic model based on equation (4) or (5) can 
be used in the aircraft equations of motion for stability 
and control studies involving either nonlinear or linear 
unsteady aerodynamics, respectively. The resulting 
equations of motion will be represented by a set of 
integro-differential equations. 

Design of Experiment 

The most important part of the experiment design is the 
selection of input forms. It has been recognized that the 
shape of an input signal could influence the accuracy of 
estimated parameters from flight measurement. Attempts 
for obtaining parameter estimates with high accuracy led 
many researchers to the development of an optimal input. 
One of the latest techniques for optimal input design is 
discussed in [6]. 

Data Compatibility Analysis 

In practice, often the measured response data, even after 
careful handling, can still contain bias errors. In order to 
verify data accuracy, a compatibility check can be applied 
to the measured aircraft responses. This check includes 
aircraft state estimation, based on known kinematics and 
the available sensor measurements, estimation of 
unknown bias errors and a comparison of reconstructed 
responses with those measured. The state equations are 
formed by kinematic relationships and the parameter 
vector usually contains constant offsets and scale factor 
errors. The estimation techniques are similar to those used 
in estimation of states and aerodynamic parameters. 

Methods of Parameter Estimation 
Model structure determination and parameter estimation 
combined with collinearity diagnostic form the principal 
part of the identification procedure. From the postulated 
model and measured data the model structure can be 
determined as explained later. When the model structure 
is known the parameter estimation can follow. In 
aeronautical application three methods, the maximum 
likelihood (ML), linear regression (LR) and extended 
Kalman filter (EKF), are the basic techniques. 

The ML estimates are obtained by maximizing the 
conditional probability of measurement 
Z=[z(l),z(2),...,z(n)] T given a value of 0, i.e. 

0 =max p(Z|0) 


Rather than minimize p(Z|0), it is more convenient to 
minimize the negative logarithm of the likelihood 
function L(0) = p(Z|0), called log-likelihood function, 

0 =min(- CnL(0)} 

0 


Substantial simplification to the ML estimation is 
obtained by assuming no external disturbances to the 
system and no measurement errors in input data. Then the 
ML estimation is reduced to the output error method with 
the cost function 

J = A Ev T (i)R' l v(i) + — fn |R| ( 6 ) 

2 i«i 2 


where v(i) are the residuals, v(i) =z(i)-y(i,0), and R is 
the measurement noise covariance matrix. Experience 
shows that a suitable technique for minimization of (6) is 
the Modified Newton-Raphson method. Using this 
technique, the step size, A0, for parameter estimates is 
given by 

A8=M~‘ (7) 

ae 


where M is the Fisher information matrix containing 
products and sums of the first order partial s of the log- 
likelihood function. The expression for the information 
matrix is 

M = (8) 

M { 00 J 00 

The information matrix provides also a lower bound on 
parameter covariances, i.e. 

Cov (0) = E{(0 - 0) (0 - 0) T } a M' 1 (9) 

The inverse of M is called the Cramer-Rao lower bound 
and the ML estimates approach it asymptotically as N 
increases. Expression (8) is valid if the measurement 
noise is random, Gaussian and white. Numerous analyses 
from measured flight data showed, however, that these 
residuals can be far from being white. 
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The residuals very often contain some deterministic 
components. The result is colored residuals leading to 
new expression for parameter covariance matrix [7] in the 
form 


Cov(0) = 




e R-« S a vv (i-j)R->M> 




j-1 


ae 


M" 


(io) 


where 9t vv (i-j) is the autocorrelation matrix for the 
output residual vector. It is estimated from 


for white noise and as 


Cov(0) = 




E x(i) E 9^(i-j)x(j) T 


j"< 


(X T X) _1 


(16) 


for colored noise where a 2 is the variance of c(i) 

^«(k) = -i- N S K e( i )e(i + k) = ^(- k ) ( 17 > 

N-K i-i 


„(k) = -i- E v(i)v(i+k) T = & vv (-k) < n > 
N-k i-i 


and 


c(i)=y(i)-x(i) T 0 


(18) 


Further simplification to the ML method is obtained by 
assuming that both the input and states variables are 
measured without errors. This assumption leads to an 
equation error method which can be formulated as a 
linear regression applied to aerodynamic model equation 
(3). When the states, and inputs in this equation are 
replaced by measured values the equivalent of the general 
regression equation 

y(i) = 0 o + 0, x i(i) + --«®n-i x »- l ( i ) +e (') (12) 

is obtained. In this equation, y now represents a 
dependent variable, Xj to x n ., are the regressors and e is 
the equation error. When the regression equations are 
expressed as 

Y = X0 e < 13 > 


the least squares (LS) parameter estimates are obtained as 

e=pt T x)-'x T Y ( 14 > 

and their covariance matrix as 

Cov(0) = o 2 (X T X)"‘ (,5) 


The linear regression is a widely used estimation method 
for the following reasons: 

a) it is a simple non-iterative method for parameter 
estimation; 

b) the LS estimates serve as nominal starting values for 
the ML and EKF methods; 

c) linear regression can be applied to data generated by 
partitioning an ensemble of data from repeated 
measurements with respect to one or more variables [8]; 

d) linear regression can be extended to a technique for 
model structure determination, e.g.stepwise regression; 

e) formulation of a regression problem can be used for 
investigation of near-linear dependence (collinearity) 
among measured state and input variables and for the 
development of biased estimation techniques for dealing 
with highly collinear data (see Section 3). 

For the development of the EKF algorithm the state 
vector is augmented by a parameter vector, x a = [x T : 0 T ] T . 
The vector x a and its covariance matrix are estimated 
from measurements by minimizing the cost function 

J=E{[x i (i|i-l)-x I (i)] T [x a (i|i-l)-x.(i)]} 

Because of inherent feed-back in the algorithm, the EKF 
can be easily applied to an unstable system for which the 
use of the output error method might be difficult or even 
impossible. On the other hand, the main disadvantage of 
the EKF method is that the initial conditions and the 
covariance matrices of process and measurement noise 
must be specified at the start of the estimation process. 
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Model Structure Determination 

A major problem in system identification is the selection, 
from measured data, of an adequate (parsimonious) 
model. An adequate model is considered to be a model 
which sufficiently fits the data, facilitates the successful 
estimation of unknown parameters whose existence can 
be substantiated, and has good prediction capabilities. For 
model structure determination a method of stepwise 
regression is being used [9]. The determination of a 
model for the aerodynamic coefficients includes three 
steps: postulation of terms that might enter the model, 
selection of an adequate model, and validation of the 
model selected. After postulating the aerodynamic model 

equations, significanttermsamongthe candidate variables 

are determined and the corresponding parameters 
estimated. At every step of the stepwise regression, the 
variables incorporated into the model in previous stages 
and a new variable entering the model are reexamined for 
their significance. Experience shows, however, that the 
model based only on the statistical significance of 
individual parameters in the model can still include too 
many terms and may have poor prediction capabilities. 
Several criteria for the selection of an adequate model 
have been, therefore, introduced. The most often used are: 

a) The computed values of F-statistic, 
given as the ratio of regression mean square to residual 
mean square. Heuristically, the model with the maximum 
F-values is the "best" one for a given set of data. 

b) The value of the coefficient of 
determination, R\ which can be interpreted as measuring 
the proportion of the variation explained by the terms 
other than 0 O in the model. 

c) The prediction sum of squares PRESS 
defined for the k* subset of model parameters as 

PRESS = 

N . . „ <' 9) 

E {y (0 - y[i | x( 1) x(i “ 1), x(i + l),.-,x(N)] k P 

i“l 

For the model to be a good predictor the value of PRESS 
should be minimal. 

Model Validation 

Model validation is the last step in the identification 
process and should be applied regardless of the 
complexity of the estimation method. The resulting model 
must demonstrate that its parameters have physically 
reasonable values and acceptable accuracy, and that the 
model is a good predictor. For those reasons the 
parameter estimates are compared with any information 
available about aircraft aerodynamics. Prediction 
capabilities of the model are checked on a set of data not 
used in the identification process. 


3. DATA COLLLINEARITY AND BIASED 
ESTIMATION 

As pointed out in the Introduction the augmentation of 
high performance aircraft very often introduces near 
linear relationships among the input and output variables 
(data collinearity). When linear regression is used in data 
analysis the collinearity results in an ill-conditioned X T X 
matrix in expression (14) for LS parameter estimates. 
Because of that the collinearity can cause computational 
problems and reduce the accuracy of estimates. Three 
procedures for detection of collinearity are recommended 
in [10] and applied to flight data. They are: 

a) examination of the correlation matrix 
X* T X* and its inverse, where the matrix X’ is formed by 
centered and scaled regressors; 

b) eigenvalue analysis of the X T X matrix 
or singular value decomposition of the X matrix; 

c) parameter variance decomposition into 
a sum of components, each corresponding to one and 
only one of the eigenvalues of the X T X matrix or singular 
values of the X matrix. 

The parameter variance decomposition approach for 
detecting collinearity was proposed in [1 ]. It follows from 

the covariance matrix of parameter estimates 0 which can 
be also obtained as 

Cov (0) = o 2 TA‘*T t (20) 

where A is a diagonal matrix whose elements are the 
eigenvalues of X T X and T is a matrix whose columns are 
the eigenvectors of X T X. 

The variance of each parameter is equal to 
2 2 

o 2 (0WE ** = o 2 S-^ (2D 

J k-1 Xj k-1 ^ 


where t k are the elements of eigenvector t associated 
with A,.. Eq. (21) decomposes the variance of each 
parameter into a sum of components, each corresponding 
to one and only one of the n singular values In (21) 
the singular values appear in denominator, so one or more 
small singular values can substantially increase the 
variance of 0 r This means that an unusually high 
proportion of the variance of two or more coefficients for 
the same small singular value can provide evidence that 
the corresponding near dependency is causing problems. 
Introducing 

t? n 

1*j ' 


the j,k variance-decomposition proportion as the 
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proportion of the variance of the jth regression coefficient 
associated with the kth components of its decomposition 
in (21) is given as 

j,k = l,2,...,n (22) 

a 


Since two or more regressors are required to create near 
dependency, then two or more variances will be adversely 
affected by high variance-decomposition proportions 
associated with a single singular value. Variance- 
decomposition proportions greater than 0.5 are 
recommended in [1] as a guidance for possible 
collinearity problems. It is also suggested that the 
columns of X should be scaled to unit length but not 
centered. Thus the role of the bias term in near-linear 
dependencies can be diagnosed. 

If the collinearity diagnostic reveals a serious problem 
some way of dealing with it should be chosen. Additional 
data can be selected, the experiment can be redesigned, 
the model can be respecified or different techniques from 
the ordinary LS procedure can be used. Several 
estimation techniques which can be applied to data with 
severe collinearity have been developed [11]. Their 
development drops the requirement that the estimator of0 
be unbiased. The new estimator, however, should have 
smaller variance than the LS estimator. By allowing the 
small amount of bias the parameter variance can be made 
small such that the mean square error of 0 is less than 
the variance of the unbiased LS estimator. Techniques 
with this property belong to a class of biased estimation 
methods. Two of the techniques, the mixed estimation 
and rank reduction regression, have been introduced in 
[10] and used in flight data analysis. The experience with 
the rank reduction regression showed some difficulties in 
its application to aerodynamic model equations, mainly 
because of small number of unknown parameters in these 
equations. For that reason, only the mixed estimation will 
be introduced. 

The mixed estimation (ME) is a procedure which uses 
prior information on parameters to augment measured 
data directly instead of through a prior distribution. 
Mixed estimation includes the usual regression model 
given by eq. (13) and the additional assumption that a set 
of prior conditions on 0 can be written as 

d = A0 + f (23) 

In this equation, A is a matrix of known constants and£ 
is a vector of random variables with E(Q-0 and 
E(i!£ T )=a 2 W, where W is a known weighing matrix. 
Combining (13) and (23) the mixed model is obtained. 


For known a 2 the application of least squares to this 
model results in the mixed estimator 

8,^ =(X t X+A t W‘A)' 1 (X T Y + A T W-'d) ( 24 ) 


It is shown in [10] that the addition of prior information 
to the ordinary regression results in reduction of 
parameter variance. In real application of the ME the a 
pri ori values are not known exactly therefore the resulting 
estimator is biased [10]. 

4. EXAMPLES 

In the following three examples the measured data from 
experiments on the X-29A, X-31A and F-16XL aircraft 
will be used. These examples will demonstrate mainly 
specific problems related to identification of high 
performance aircraft. In addition, variations of some 
parameters with the angle of attack or Mach number, and 
their correlation with wind tunnel data and flight results 
from scale model will be also shown. 

4.1 X-29A Aircraft 

The test vehicle is a single engine, single seat fighter-type 
research aircraft with forward-swept wings. The aircraft 
has highly relaxed static longitudinal stability in subsonic 
and transonic regimes and near-neutral stability in 
supersonic regimes. For longitudinal control, deflections 
of canard, wing flap (flaperon), and fuselage strake are 
used. The lateral control is provided by the rudder and 
asymmetric deflection of flaperon. In addition to manual 
control of the aircraft, the concept of remotely augmented 
vehicle (RAV) could be used for the excitation of aircraft 
responses. The RAV arrangement employs a ground 
computer to augment the onboard control system. This 
capability is used to introduce a command to the control 
stick (pitch stick or roll stick command), rudder pedal, or 
individual control surfaces. The RAV commands, usually 
a pulse or doublet, are summed onto the already existing 
commands in order to independently move flaps, strake, 
canard, rudder, or differential flap. More about the RAV 
system can be found in [12]. A drawing of the aircraft is 
presented in figure 2. A more detailed description of the 
aircraft and its control system is contained in [13]. 

The following three sets of data were available for 
estimation of aircraft parameters: 

1. longitudinal maneuvers excited by a pilot at Mach 
numbers from 0.5 to 1.4, 

2. longitudinal and lateral maneuvers with computer 
generated inputs (RAV experiment) at Mach numbers 
from 0.6 to 1.3, 

3. low speed lateral maneuvers initiated by a pilot at the 
angles of attack between 8° to 50°. 
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During data analysis several problems associated with 
inherent instability, high augmentation and sometimes 
insufficient excitation of maneuvers had to be addressed. 
Among them were: 

1. parameter estimation of an unstable vehicle, 

2. data collinearity and its diagnostics, 

3. adverse effect of data collinearity on parameter 
identifiability and accuracy. 

Because of these problems, a linear regression and mixed 
estimation were used in data analysis. 

Time histories of a typical longitudinal pilot and 
computer-generated input are presented in figures 3 and 

4. From figure 3 close correlation among all open-loop 
inputs is obvious. The change in data collinearity caused 
by replacing the pilot by computer-generated input as a 
sequence of commanded flap, strake, canard and stick 
deflections is demonstrated in figure 4. The collinearity 
was assessed by comparing correlations between 
regressors, and corresponding, condition indexes and 
parameter variance proportions. The condition index is 
defined as the ratio of the maximum eigenvalue of the 
information matrix X T X to one of the remaining 
eigenvalues. 

The correlation matrices for the two sets of data are given 
in tables I and II, respectively. The data with the pilot 
input show correlation between 
5 c and qc/2Vand 5 c and 5 s . The RAV data do not exhibit 
any significant correlation between regressors. The 
condition indexes and parameter variance proportions are 
given in tables III and IV. In this example the maximum 
condition index (condition number) for the set of data 
with pilot input is 174, in the second case 14, thus 
indicating reduced spread of eigenvalues where the RAV 
system was used. The variance proportions in table III for 
the largest condition index show r strong collinearity 
among the bias term, canard and strake effectiveness. The 
same quantities in table IV indicate only a possibility of 
collinearity between the bias term and canard 
effectiveness. 

Table V demonstrates a possibility of estimating 
parameters in the regression equation for the pitching- 
moment coefficient with sufficient accuracy. Included are 
the increments in coefficient of determination, AR 2 , and 

t-statistics, t*. The values of AR 2 represent the amount of 
information in the data explained by the individual terms 
in the model, t-statistic can be considered as a measure of 
significance of individual parameters. The data with the 
pilot input revealed that C mfic 5 c is the highly influential 
term in the pitching-moment equation, that there is a 
limited possibility for accurate estimates of parameters 

C ro5f and C^, and that the significance of theC^qc/W 
term is almost zero. The RAV experiment improved the 
identifiability of parameters C ma . C m6f and C m6s , still 
maintaining the C m6c S c as a dominant term. The chance 


for accurate estimation of C mq remains small. 

The estimates of two parameters C ma and C m6c which 
contribute the most to the pitching moment are plotted 
against the Mach number in figures 5 and 6 respectively 
and compared with the wind tunnel results. The estimates 
were obtained from data with pilot input using the least 
squares and mixed estimation, and from RAV experiment 
using only the least squares method. The wind tunnel 
values of the strake effectiveness were used as a priori 
values. The accuracies of the a priori values were 
determined from repeated test in two different wind 
tunnels. As can be seen from these figures, the accuracy 
of parameters was improved either by applying the mixed 
estimation to data with pilot input or by using data from 
the RAV experiment. Similar conclusion can be drawn 
from average standard errors in estimated parameters. 
These values are shown in table VI together with the fit 
error for the coefficient C m . 

The measured data from all low speed lateral maneuvers 
were assembled into one set with 51,200 data points 
which was then partitioned into 42 one-degree- a subsets 
and 1 three-degree- a subset. A half of the selected lateral 
maneuvers was analyzed as individual maneuvers using 
stepwise regression. The possibility of data collinearity in 
measured data was investigated by procedures explained 
in the previous examples. Application of stepwise 
regression to partitioned data resulted in models for the 
lateral aerodynamic coefficients and least squares 
estimates of parameters in these models. For the data 
subsets with a<40°, models with linear stability and 

control derivatives were adequate. For data at a>40 () . 
some nonlinear and longitudinal control terms were 
selected by the stepwise regression. These additional 
terms did not provide any comprehensive information 
about aerodynamic nonlinearities or effects of 
longitudinal control setting on the lateral aerodynamic 
coefficients. 

The estimated parameters (stability and control 
derivatives) from 43 subsets were plotted against the 
angle of attack and fitted by quadratic polynomial splines. 
In addition to fitted splines, the 2 ct - confidence limits on 
the mean were computed as 


0(a)±2s v / x T (X T X) 1 x 


In (25), 0 is the mean value of a parameter given by the 

fitted spline, s is the standard error of 0 estimated from 
the residuals, x is a vector of regressors. As an example, 
estimated values of the parameter C np , fitted spline, and 
2 ct - confidence limits are shown in figure 7. 

Stepwise regression analysis of data from single 
maneuvers showed that linear models for the aerodynamic 
coefficients were adequate within the angle-of attack 
range from 8° to 40°. For angles of attack greater than 40° 
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models for the lateral-force and yawing-moment 
coefficient included some of the nonlinear terms 

p 2 , p 3 , aP, p6 c or 8 r 3 . As for the partitioned data, there 
was no consistent information about parameters associated 
with these nonlinear terms. 

The least squares parameter estimates at low and 
moderate angles of attack were close to results from 
partitioned data. For a>20°, however, increased scatter 
in the estimates and their deviation from the previous 
results using partitioned data were observed. This 
inconsistency was caused by data collinearity detected 
mostly among the variables pb/2V, 5 a and 8 r . In applying 
the mixed estimation technique, the a priori values were 
selected for parameters which were affected by 
collinearity. Their changes had only a small effect on 
aerodynamic coefficients. The a priori values and their 
uncertainty were set to the mean values and their standard 
errors obtained from partitioned data. An example of data 
collinearity is given in table VII where the eigervalues of 
the information matrix, condition indexes and parameter 
variance decompositions are included. The damaging 
effect of data collinearity on parameters associated with 
the variables pb/2V, 8 t and 8 r can be expected. Table 
VIII presents the parameter estimates and their standard 
errors in the rolling-moment equation using stepwise 
regression and, mixed estimation. In the last column of 
table VIII the parameter estimates from partitioned data 
are given. Large changes in the least square estimates of 
C, p , C r6a and C f6r are visible during the model 
development. The final estimates of the three parameters 
are also quite different from those obtained from 
partitioned data. On the other hand, the parameters from 
the mixed estimation technique are close to the result 
from partitioned data. These estimates also have lower 
standard errors than the least squares results. The decrease 
in the fit error, s (C f ) , and squared correlation coefficient, 
R 2 , are not substantial (see last two rows in table VIII). 

Selected parameters obtained from partitined data and 
single maneuvers are plotted against the angle of attack 
in figures 8 and 9, and compared with wind tunnel data. 
The 2a - confidence limits were omitted in these figures. 
However, the minimum and maximum values for standard 
errors of the estimated parameters from partitioned data 
are given in table IX. Both figures indicate no significant 
differences between the two sets of flight results. Large 
differences, however exist between flight and wind tunnel 
results in figure 8. Flight data exhibit a sudden increase 
in parameter values of C >p and C np at angles of attack 
around 40° to 45°. As indicated by wind tunnel 
investigations in [14], an increase of lateral force and 
yawing moment due to sideslip is caused by the forebody 
vortex assymetry. This assymetry can produce a sideforce 
which moves the nose into the sideslip and thus enhances 
directional stability. The same effect was observed in [15] 
during wind tunnel testing of the scaled model. The 
present differences between flight and wind tunnel results 
can be caused by different Reynolds numbers in the two 
experiements (0.68 x 10 6 in the wind tunnel; approx. 6 x 
10 6 during the flight test) and by the effect of the nose 


boom on the full scale aircraft Positive yaw damping for a >30° 
is predicted by flight data whereas the wind tunnel results 
show low yaw damping over the whole range of angles 
of attack. In figure 9 three rolling moment parameters are 
included. The parameter C rp estimated from flight data 
agrees, in general, with wind tunnel predicitions. The roll 
damping decreases above a = 15°, from wind tunnel 
above a =25°. Positive values of C lp from flight are 
about three times higher than those from wind tunnel. As 
mentioned in [15] and [16] the forebody aerodynamics 
dominates the roll-damping parameter at high angles of 
attack and causes the unstable damping. Therefore the 
differences between flight and wind tunnel results maybe 
attributed to different forebody aerodynamics in these two 
test conditions. The parameter C /6a is estimated with high 
consistency and agrees well with wind tunnel data for the 
angles of attack between 12° and 40°. Some differences 
exists outside this interval. More about the parameters can 
be found in [17]. 

The results from this example lead to the following 
conclusions: 

1. A collinearity diagnostic and a comparison of 
parameter estimates from flight data using the RAV 
experiment with those from flight data using pilot input 
revealed that the computer-generated deflections of 
individual control surfaces can substantially decrease data 
collinearity. A simple least squares technique can be used 
in estimation of all parameters in the model which also 
means estimation of the effectiveness of all controls used. 

2. The experiment providing the data for the analysis was 
not properly designed. The selected input forms resulted 
in very small excitation of response variables. As a result 
some parameters were not identifiable and the 
identifiability of several remaining parameters was 
significantly reduced. These identifiability problems were 
apparent from the diagnostic of selected maneuvers. 

3. For low speed lateral maneuvers linear aerodynamic 
models determined by applying stepwise regression 
techniques to partitioned data and single maneuvers were 
found to be adequate for the angles of attack less than 
40°. For angles of attack greater than 40°, nonlinear and 
longitudinal control terms entered the models selected by 
the estimation technique. These terms, however, did not 
provide any comprehensive information about 
aerodynamic nonlinearities and the effect of longitudinal 
control on lateral aerodynamic coefficients. Because of 
data collinearity detected in single maneuvers, the data 
from these maneuvers were reanalyzed by using mixed 
estimation. No significant differences existed between 
parameters estimated from partitioned data and those 
estimated from single maneuvers. 

4.2. X-31A Aircraft 

The X-31A is a single engine, single seat fighter-type 
research aircraft with delta wing and canard. The aircraft 
is inherently unstable in subsonic regimes and is 
controllable by a digital control system. For longitudinal 
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control, canard and symmetrical flaperon deflections are 
used. The lateral control is provided by the rudder and 
differential flaperon. In addition to the aerodynamic 
control surfaces, a three-vane thrust vectoring system is 
mounted around the engine exhaust nozzle. This system 
allows thrust deflection up to 15° and is used for 
augmentation of pitch and yaw control during low speed 
and post-stall flights. A drawing of the aircraft is 
presented in figure 10. 

Maneuvers from eight flights were selected for 
aerodynamic model structure determination and parameter 
estimation. These maneuvers were initiated from trim 
conditions at altitudes between 6,000 to 9,000 meters and 
angles of attack between 10° and 70°. Two types of input 
were used, the pilot inputs and inputs generated by a 
flutter test box installed on one of the aircraft. The pilot 
inputs were either pitch command or yaw and roll 
commands in the form of a single doublet or a 
combination of doublets. The flutter test box allowed 
separate excitation of all aerodynamic control surfaces. 
These inputs were in the form of ”3211 ” multiple signals. 
The surface deflection limits in the tests were 7. 1° for the 
canard, 11.8° for the flaperons and 7.8° for the rudder. To 
ensure good excitation of aircraft motion, the longitudinal 
inputs were in some cases preceded by the pilot pitch 
doublets and the lateral inputs by pilot roll and yaw 
doublets. Maneuvers with pilot inputs were flown with 

thrust vectoring off (a<30°) and on. The single surface 
excitation maneuvers were realized only with thrust 
vectoring on. 

Estimated parameters from flight data were correlated 
with those obtained from wind tunnel data and drop 
model experiment. Static and dynamic wind-tunnel test 
were conducted at the 30 - by 60 - Foot Tunnel at NASA 
LaRC using a 1 9 - percent scale model (see [ 1 8] for some 
results). For a comparison with flight results presented, 
the parameters C fp and C fr were estimated from the 
oscillatory data using the techniques of [19] and [3]. 
More about this analysis will be presented in the third 
example. 

In preparation for X-31A flight testing, the unpowered, 
27-percent dynamically-scaled radio-controlled model of 
the X-31A was built and tested at NASA LaRC. During 
the test, the model was attached to a helicopter and lifted 
to an altitude between 1,500 and 3,000 m, and then 
released with the helicopter in forward motion. At an 
altitude of about 300 m the flight was terminated by 
deploying a large parachute. The model was controlled by 
a pilot on the ground but the model also had the 
capability to accept preprogrammed surface deflection 
commands. This feature was used in generating 
maneuvers for the purpose of parameter estimation. 
During the experiment, only the lateral parameters were 
estimated. They are presented in [20]. In this paper only 
estimates of four lateral parameters and two thrust 
vectoring parameters are presented. More results can be 
found in [21]. 


Lateral parameters: 

In maneuvers with pilot inputs, high augmentation 
prevented sufficient excitation of sideslipangle and lateral 

acceleration, especially at angles of attack above 30°, and 
introduced high pairwise correlation between variables p 
and r, 5 a and 8 r , 5 r and 6 qtv . The effect of limited 
excitation and high correlation was reflected in low 
accuracy of parameter estimates and, in many cases, 
inability to estimate some parameters at all. In the 
rolling-moment equation the terms with 5 a and P were 
predominant, explaining about 75% variation in the data 
of the 87% for the complete model. The main influential 
terms in the yawing-moment equation were those with 
either 8 r (for low a) or 5 qtv in maneuvers at high angles 
of attack. Their contribution represented about 85 to 96% 
from the overall value of 98% for the model with all 
terms included. 


The parameter identifiability was improved, in general, by 
the introduction of a single surface excitation 
(5 t or 8 r inputs) . The above mentioned correlations were 
substantially reduced and the amplitude of sideslip angle 
increased. Asa result of that, parameter identifiability in 
the rolling-moment equation was improved. However, the 
problem of accurate parameter estimation in the yawing- 
moment and lateral-force equations remained, with the 
exception of the thrust-vectoring effectiveness. Four 
parameters in the rolling-moment equation 


= 

<VC, p P + C, p |J + C,4 + C (St 6, + C (6t 6 t 


(26) 


are presented in figure 1 1 with their 2a -error bounds. In 
this figure the results from wind tunnel and drop model 
tests are also included. The estimates of C rp are 
consistent with the exception of the region around 
a =25°. There is significant departure of flight results 
from wind tunnel data for the angle of attack between 30 
and 45°. The reason for that has not been explained yet. 
Small scatter is also apparent in the estimates of C ?&a 
The flight data demonstrate about 50% reduction in 
aileron effectiveness when compared to wind tunnel 
values. The estimates of C rp and C, r demonstrate the 
advantage of single surface excitation over pilot induced 
maneuvers. In the latter case, high correlation between 
rolling and yawing velocities, and low information 
content in the data prevented the estimates of both 
oscillatory parameters. There are no substantial 
differences between full-scale flight and wind tunnel 
estimates of these two dynamic parameters. The drop 

model parameter values at 30°<a <45° are scattered with 
values lower than indicated by results from full-scale 
aircraft test. 

Thrust vectoring effectiveness: 

The parameter C m6qtv was estimated directly from 
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maneuvers with the 8 c inputs because the adequate model 
for the pitching moment included only variable 
a,6 e and 6^. These three regressors were not 
correlated. The maneuvers with 6 c inputs provided 

estimates of for a <45° and C^ 5qtv for a >45°. 
From these estimates the thrust vectoring effectiveness 
was computed as 


or 


“ni&qtv 


=( c '» 



c 

mfiqtv mftqtv 



(27) 


where the relation 8 e / 8 qtv was evaluated from measured 
data and where the values of C mftc were obtained from 
maneuvers with pilot input (no thrust vectoring) and5 c 
input. The values of the thrust-vectoring effectiveness, 
Cn 5 itv* were estimated from the data with pilot input and 
single surface excitation. In the first cases, only 
maneuvers for a >35° were used because the effect of 
rudder deflection on the yawing moment was small. 
When necessary the final estimates were corrected for the 
effect of the rudder deflection by computing C n8rtv from 
the estimated values of C " n6ftv as 


"'nfirtv 


= c \ 



(28) 


From maneuvers excited by a single surface deflection, 
the parameter C n6rtv was estimated directly. Both 
parameters C m5qtv and C n6rtv are plotted against the thrust 
coefficient in figure 12 and compared with their 
theoretical values 


C 


mfiqtv 




(29) 


where x ft is the distance of thrust impact point from 
airplane center of gravity. So far there has been no 
explanation for the differences between the theory and the 
experiment. 

From the results presented here and in [22] the following 
conclusions can be drawn: 

1. In piloted maneuvers, high correlations between input 
and output variables were observed. These correlations 
resulted in low accuracy of estimated parameters and 
prevented estimation of effectiveness of several control 
surfaces. 


parameters. 

3. The accuracy of parameters was also affected by their 
low sensitivity and small excitation of response variables, 
especially at high angles of attack. 

4. All the influential parameters were, in general, in 
agreement with wind tunnel data and results from drop 
model tests. Some unexplained differences were, 
however, observed. 

5. The estimated thrust vectoring effectiveness was found 
lower than its theoretical values. 

6. The predictive capabilities of the resulting model 
determined from flight test data are very good for low 
amplitude maneuver at low to moderate angles of attack. 
Some deterioration in prediction was observed at high 
angles of attack. 

4.3. F-I6XL Aircraft 

For better understanding of aircraft aerodynamics in large 
amplitude maneuvers, NASA LaRC conducted a series of 
wind tunnel tests on a model of the F-I6XL. A sketch of 
the 0. 1 0 - scale model is shown in figure 1 3 together with 
some basic model dimensions. Tests included 
measurement of aerodynamic forces and moments under 
static conditions, followed by oscillatory test about all 
three body axes and ramp tests in pitch. One of the 
reasons for the testing was to determine a mathematical 
model with unsteady aerodynamic terms from oscillatory 
data at different angles of attack and frequencies. In this 
paper only limited results from small amplitude (±5°) 
oscillations in pitch will be shown. More results can be 
found in [3] and upcoming reports. 

For the following analysis of the oscillatory data it is 
assumed that the aerodynamic coefficients are linear 
functions of the angle of attack, pitching velocity and 
their rates. Then, for example, the increment in the lift 
coefficient with respect to its mean value can be 
formulated as 

ac l = 

Q B ( C\ 2 ( 30 ) 

c^q 


then for the harmonic motion 

Aq = «A ( C L.' k 2 C L 4 ) SinG,t + a A k ( C Li +C L,) COS “ t 

=a A (C L .sin w t + kCj^cos u> t) 

(31) 

where a A is the amplitude, © the angular frequency and 
a AC L . =a A( C L.- kJC L4) 


2. Single surface excitation reduced the above mentioned 
correlations. It was possible to estimate the effectiveness 
of all controls and increase the accuracy of estimated 


«A kC L, = a A k ( C L i +C U 1 ) 
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represent the Fourier coefficients. The in-phase and out- 
of-phase components of C L (C L and C L ) can be obtained 

by integrating the time histories of AC L over a selected 
number of cycles. The in-phase and out-of-phase 
components of C L are plotted against the angle of attack 
in figure 14. The figure shows the effect of frequency 
which is especially strong on the out-of-phase 
components. A comparison of steady and oscillatory data 
is given in figure 15. 

A strong dependence of wind tunnel oscillatory data on 
frequency lead to the development of models with 
unsteady aerodynamic terms, see e.g. [20] and [3]. As an 
example, the model for the lift increment was formulated 
in [20] as 

ACl= 

(32> 

0 0 

where and <2^(0 are the indicial functions. 

For obtaining the model with limited number of 
parameters, it was assumed that the effect of q(t) on the 
lift could be neglected and the indicial function C u (t) 
can be expressed as 

C La (t)=a(l-e' b, Vc ( 33 > 

Considering the above mentioned assumptions, equation 
(33) is simplified as 

ac l = 

* 

C^Waft) -aj‘e' b,<, ' ,) -^-a(T)dT +^C Lq («)q(t) 

(34) 

where C La (<») and are the rates of change of 

C L with a and q in steady flow. The steady form of 
equation (32) for harmonic changes in a is identical to 
that of equation (31), that is 

AC L (t) = a A C Lo sina)t + a A kC Lq coso)t (^5) 

However, with the indicial function of equation (33), the 
expressions for and have the form 

x 2 k 2 

(36) 

1 +T 2 k 2 


From the experiment, the in-phase and out-of phase 
components are usually obtained for different values of 
the angle of attack and reduced frequency while keeping 
the amplitude of the oscillations constant. Then the 
equations (36) and (37) can be generalized as 


“ji =u r a i z uj 

(38) 

V v i _a i 2 vi 

(39) 

These equation define Model I where 
coefficient 

for the lift 

U i =C L«( a i) V i =C Lq<“i) 


. _ T i*S . . *. 


“ J 1+T^k* VJ 


i = l,2,...,n 
for 

j = 1,2,... an 


In equations (38) and (39) there are, in general. 3n + 1 
unknown parameters: u^v^a and t,. They can be 
estimated from experimental data u j; and v.. by 
minimizing the cost function 

J ,= E s([v(V¥ U) )fHv(vr 

j-1 1-1 

-Vvj)] 2 } 

In formulating airplane equations of motion it might be 
more convenient to obtain expressions for u, v, and a as 
a function of the angle of attack rather than their discrete 
values. For that reason the previous model was 
reformulated as Model II defined by the following 
equations: 

S ji =u(a i )-a(a i )z uj 

(40) 

Vj, = v (<*j) ~a(“i)z V j 

(41) 


The form of expressions for u(a) y v(a) and a(a s ) can 
be either specified from the variation of estimated 
parameters in Model I or the form can be postulated in 
terms of polynomials and/or polynomial splines. In the 
second case adequate models for u(a),v(a) and a(a) can 
be determined from measured data by a stepwise 
regression. The previously estimated value of parameter x, 
can be used as an a priori value thus simplifying the 
parameter estimation procedure using Model II and the 
cost function 


1 +T,k ii 


(37) 


where x, = V/b, i is the nondimensional time constant. 


J n =2 

j-1 i-1 

+[v ji -(v(« i )-a(a j )z vj )] 2 } 


Both estimation procedures were applied to oscillatory 
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data at four frequencies. The data for frequency k = 0. 19 
were used to demonstrate model prediction capabilities. 
The parameters C La (~), C^O*) and a in both models are 
plotted in figure 16. The agreement between both sets of 
results is very good. The parameter C La («) agree well 
with that from static data. Low accuracy in C^*) could 
result from small number of data points and small 
sensitivity of this parameter. The parameter a indicated 
smooth variation of unsteady term with the angle of 
attack and the largest effect of unsteady terms on the 
coefficient at a around 40°. The estimated value oftj 
was 17.2±1.0, which means b, =2.71 ±0.16(sec). This 
result indicates that the time constant associated with the 
unsteady effect is about 0.4 sec. The predicted 
components from model eq.(40) and (41) for k =0.1 9 are 
presented in figure 17 together with the corresponding 
measured values. This figure demonstrates that Model II 
is a good predictor. Similar result was obtained for Model 
I. 

As follows from this example, a strong dependence of 
wind tunnel oscillatory data on frequency led to the 
development of model with unsteady aerodynamic terms 
in the form of indicial functions. These functions were 
postulated as simple exponentials where the unknown 
parameters included aerodynamic derivatives, the 
exponent and multiplication term. 

5. CONCLUDING REMARKS 

System identification applied to aircraft has proven to be 
a powerful tool for aircraft modeling based on 
experimental data. In recent years the introduction of 
highly maneuverable, inherently unstable and highly 
augmented aircraft has required high sophistication in 
aircraft identification methodology. At present this 
methodology includes basic steps: model postulation, 
experiment design, data compatibility analysis, model 
structure determination and parameter estimation 
combined with coll inearity and identifiability diagnostics, 
and model validation. 

The examples presented demonstrated the need for the 
introduction of an on-board computer generated input 
signals. This system and the optimal input design 
procedure improve identifiability of aerodynamic models 
and accuracy of estimated parameters. The various 
modifications of the ordinary least square method proved 
to be basic tools for model structure determination, 
parameter estimation and data collinearity assessment. 
The resulting parameters can be used as nominal values 
in the application of two parameter estimation techniques 
used in aircraft application, the maximum likelihood 
method (usually in the form of output error method) and 
the extended Kalman filter method. For meaningful 
information about the identified model, the parameter 
estimates should be completed by their covariances thus 
indicating the accuracy of parameters and their 
correlation. 


conditions were available. By repeating the identification 
using data at different trim conditions, an extended model 
could be constructed. However, for obtaining a global 
model from flight data a dedicated experiment would be 
needed. In this experiment modeling and optimal input 
design would play the major role. One of the possibilities 
for an extension of traditional aerodynamic model is the 
inclusion of unsteady effects. Some initial work has 
already been done mostly using data from dynamic wind 
tunnel testing. A strong dependency of oscillatory data on 
frequency led to identified models with time varying 
terms in the form of indicial functions. Based on these 
experiences the extension to flight data will follow. 
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Table I. X T X matrix in correlation 
form. Pilot input. 


a 

qc_ 

2V 

5, 

5, 

5c 

1.000 

0.151 

0.753 

0.711 

mszm 


1.000 

0.318 


-0.979 1 



1.000 

0.643 





1.000 

■»*£!■ 





MEM 


Table II. X T X matrix in correlation 
form. RAV experiment. 


a 

qc_ 

2V 

5, 

5, 

5 C 

1.000 


0.673 

0.328 

gtlriritM 


1.000 

0.345 

-0.375 

-0.089 



1. 000 

0.098 

-0.578 




1.000 

-0.303 





1.000 
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Table HI. Collinearity diagnostic. Pilot input. 


Eigenvalue 

Condition 

index 

Variance proportions (scaled regressors) 

l 

a 

SE- 

iv 

5, 

Sr 

5 C 

3.3102 

1 

0.0000 

0.0092 

0.4806 

0.0017 

0.0006 

0.000 i 

1.2781 

3 

0.0099 

0.000 1 

0.0397 

0.0720 

0.0058 

0.0001 

1.0333 

3 

0.0001 

0.2572 

0.1604 

0.0000 

0.0342 

0.0002 

0.2613 

13 

0.0504 

0.1140 

0.0991 

0.0762 

0.1887 1 

0.0245 

0.0981 

34 

0.1220 

0.6194 

0.2143 

0.0041 

0.7123 

0.0548 

1 0.0190 

174 

0.8176 

0.0001 

0.0061 

0.8460 

0.0585 

0.9204 


Table IV. Collinearity diagnostic for RAV experiment. 


Eigenvalue 

Condition 

index 

Variance proportions (scaled regressors) I 

i 

a 

SE- 

TS 

6, 

Sr 

Sc 

2.5320 

i 

0.0003 

0.0017 

0.8594 

0.0003 

0.0000 

0.0018 

1.4339 

2 

0.0905 

0.0002 

0.0049 

0.0008 

0.0132 

0.1901 

1.0028 

3 

0.0156 

0.2543 

0.0098 

0.1481 

0.1685 

0.2565 

0.4745 

5 

0.2173 

0.0355 

0.0171 

0.0178 

0.4886 

0.0172 

0.3798 

7 

0.0460 

0.6463 

0.0220 

0.4605 

0.0094 

0.0506 

0.1771 

14 

0.6302 

0.0621 

0.0867 

0.3724 

0.4720 

0.7147 


Table VI. Average standard errors of 

Table V. Identifiability diagnostic. estimated parameters. 


Parameter 

RAV 

Pilot input j 

LS 

LS 

ME 


0.029 

(0.014) 

0.058 

(0.075) 

0.054 

(0.020) 

S«^) 

— 

0.090 

(0.037) 

0.048 

(0.020) 

s(C M ) 

m 

■H 



Note: figures in parentheses are standard errors. 


Parameter 

RAV 

Pilot input | 

AR 1 . % 

lt*l 

WZSEM 

it*i 


9.09 


5.06 

28.0 

■9 

1.08 

13.7 

0.01 

0.0 

C -W 

75.74 

99.1 

91.49 

12.3 

■91 

18.35 

46.4 

0.94 

9.2 

LS* 

4.46 

42.7 

1.00 

9.2 


Table VII. Collinearity diagnostic. 


Eigenvalue 

Condition 

index 

Variance proportions (scaled and centered regressors) j 

p 

pb_ 

TS 

rb 

TS 

5» 

Sr 

2.387 

1.0 

0.000 

0.005 

0.002 

0.013 

0.007 

1.849 

1.3 

0.085 

0.000 

0.084 

0.006 

0.004 

0.537 

4.4 

0.167 

0.001 

0.528 

0.022 

0.001 

0.209 

11.4 

0.540 

0.000 

0.002 

0.207 

0.079 

0.017 

140.4 

0.208 

0.994 

0.384 

0.752 

0.9t0 
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Table VIII. Identifiability diagnostic and parameter estimates. 


Parameter 

Stei 

pwise Regression 

Mixed* 

estimation 

Partitioned 

data 

n = 2 

n = 3 

n = 4 

n = 5 

n = 6 

C, 

-0.0036 

-0.0014 

-0.0028 

-0.0031 

-0.0020 

-0.0172 

mbSm 

c ', 

■HI 

■HI 


■ 

IBi 

m 

HI 


- 

- 

■n 

mm 

mm 

m 


C K 

- 

- 

- 

- 

n 


■ 


- 

0.028 

(0.0058) 

Hi 

m 

| 

0.106 

(0.0046) 

0.064 

(0.0035) 

0.064 

(0.0041) 


- 


- 

0.031 

(0.0053) 


0.015 

(0.0010) 

lrT*4ftl 

mzssm 

Wm 

0.0045 

86.8 

0.0043 

87.6 

0.0040 

89.7 

■ 

■ I 

0.0040 

88.0 

- 


* a priori value C,^ = 0.010 ± 0.0021 


Note: figures in parentheses are standard errors 


Table IX. Standard errors of estimated 
parameters from partitioned data. 


Parameter 

Standard Error 1 

min 

max 


0.045 

0.15 

C r, 

0.42 

1.1 

C Y, 

0.78 

2.3 

Cy> 

0.026 

0,060 

C r. 

0.016 

0.053 

c ‘, 

0.0076 

0.016 

c >. 

0.040 

0.11 

c. 

0.063 

0.18 


0.0037 

0.01 1 

C » 

0.0021 

0.0044 

C ., 

0.0064 

0.020 


0.054 

0.12 

Q 

0.086 

0.31 


0.0029 

0.0070 


0.0013 

0.0051 



Figure 1. Block diagram of aircraft identification. 



Figure 2. Drawing of X-29A aircraft. 
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Figure 3. Time histories of measured input variables. 
Pilot input. 
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Figure 4. Time histories of measured input variables. 
RAV experiment. 



Figyre 5. Static-stability parametr estimated from 
flight data with pilot input by using (a) least 
squares, (b) mixed estimation, (c) from RAV 
experiment. 
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